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Using a sum rule approach we calculate the frequency shifts 
of the quadrupole oscillations of a harmonically trapped Bose 
gas due to the presence of a quantized vortex. Analytic re- 
sults are obtained for positive scattering lengths and large 
A'^ where the shift relative to excitations of opposite angular 
momentum is found to be proportional to the quantum cir- 
culation of the vortex and to decrease as N~^^^. Results are 
also given for smaller values of A'^ covering the transition be- 
tween the ideal gas and the Thomas-Fermi limit. For negative 
scattering lengths we predict a macroscopic instability of the 
vortex. The splitting of the collective frequencies in toroidal 
configurations is also discussed. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80. Pj, 67.40.Db 



After the experimental realization of Bose-Einstein con- 
densation in dilute atomic gases, the study of the coUe- 
cive excitations of these unique inhomogeneus quantum 
systems has been the object of both experimental 0-^] 
and theoretical work (for an update list of references 
see, for example, [0). These oscillations are characterized 
by proper quantum numbers, reflecting the symmetry of 
the confining potential. In an axially symmetric trap 
the third component of angular momentum is a natural 
quantum number and if the system is in a time reversal 
invariant configuration, elementary excitations carrying 
opposite angular momentum are degenerate. This de- 
generacy is in general removed if time reversal symmetry 
is broken. The purpose of this work is to describe the 
frequency shifts produced by the presence of a quantized 
vortex. In view of the important role played by vortices in 
understanding the mechanisms of superfluidity, the pos- 
sibility of their spectroscopic diagnostics is highly inter- 
esting since the measurements of collective frequencies 
can be carried out with high precision in these systems. 

The occurrence of splitting in the presence of a vor- 
tex can be simply understood by noting that the aver- 
age velocity flow associated with the collective oscillation 
can be either parallel or opposite to the vortex flow, de- 
pending on the sign of the angular momentum carried 
by the excitation. This produces a shift of the collective 
frequency of order Su/lo v/c where v 1/R is the 
velocity of the vortex flow and c is the sound velocity. In 
a trapped Bose gas c increases linearly with the radius 
R of the condensate, while lu is practically independent 
of R, so one expects relative shifts of the order of 1/R? . 



These effects are larger than the typical corrections to 
the Thomas-Fermi limit due to finite size effects which, 
in the absence of the vortex, behave like logR/R^ 
The problem of the frequency shift produced by a quan- 
tized vortex has been already the object of theoretical 
work using semiclassical approaches based on a large N 
expansion , as well as by full numerical solution of the 
linearized equations of motion px| ]. 

In this letter we develop a sum rule approach which 
is expected to provide exact results for the splitting of the 
excitation frequencies in large systems. The method can 
be also applied to calculate the frequency shifts for small 
values of N as well as for negative scattering lengths. Let 
us introduce the strength distribution function 



S±{E) = ^ |(n |F±| Q)\H{E - huj^o) 



(1) 



relative to the operators F± — '^'^=1 f^i^k) carrymg 
opposite angular momenta. In equation (M hujno — 
(En — Eo) are the excitation energies relative to the eigen- 
states In) of the Hamiltonian 
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which describes N interacting bosons confined by an ex- 
ternal potential. The potential Vcxt{^) — Miuj'jj'^ + 
tj^z^)/2, with r\ = + , is assumed to be axially sym- 
metric, and the interatomic force is a contact 2-body in- 
teraction whose coupling constant g = Anh^a/M is fixed 
by the s-wave scattering length a. 

In the following we will focus on the collective oscilla- 
tions of low multipolarity which are easily excited in ex- 
periments by suitable modulation of the harmonic trap. 
For the quadrupole case we will consider the modes ex- 
cited by the operators 



and 



/± = (a; ± iyf 



/± = (x ± iy)z 



(3) 



(4) 



carrying angular momentum m = ±2 and m = ±1 re- 
spectively (here and in the following we will identify m 
with the third component of angular momentum of the 
elementary excitation). Only excitations with m ^ are 
relevant for the present discussion. 
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In the absence of vortices the ground state has zero an- 
gular momentum and, for large N and positive scattering 
lengths, the collective states excited by the operators (||) 
and are well described by hydrodynamic theory of 
superfluids. This yields the result uj± = y/2uj± and 
Lu±= -y/ oj'^ + ljjI for the m = ±2 and m = ±l frequencies. 
Notice that these results differ from the ideal gas predic- 
tions, tjO±=2ijj±_ and U!±=u!±_ -\- uiz, as a consequence of 
interaction effects which suppress the contribution of the 
kinetic energy pressure term in the equations of motion. 
Differently from the quadrupole excitations, the frequen- 
cies of the dipole modes excited by f± — x ± iy are in- 
stead unaffected by two body interactions and are given 
by w± — uj±. This behavior is the consequence of the 
translational invariance of the interatomic force which 
cannot affect the motion of the center of mass, even in 
the presence of a vortex. 

The moments 



the sum rule vanishes. Notice that results (|1 
are formally independent of the choice of the external 
potential as well as of the two-body interaction. 



The results for the moments 



and 



be 
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used to calculate the shift of the collective frequencies 
when the number of atoms in the trap is large and a 
is positive. In fact in this limit, where the behavior of 
the system is properly described by hydrodynamic theory 
of superfluids, one expects that the moments calculated 
above will be exhausted by two modes with frequency uj± 
excited, respectively, by the operators F± : 



(12) 



m^= dE{S+{E)±S-{E))EP 
Jo 



(5) 



of the strength distribution can be calculated using 
closure relations. For the lowest moments we find the 
results: 



2^ = {[F^,F+])=0 

.+ = ([F_,[i?,F+]]) = — (|VA 



(6) 
(7) 



where cr^ are the corresponding strengths. Assumption 
( |l^ ) is equivalent to a Bijl-Feynman ansatz often 
used to describe the collective excitations in interacting 
many-body systems. In the case of superfluid helium it 
provides an exact description of the excitation spectrum 
in the phonon regime (for a recent discussion on sum 
rules and collective excitations in Bose superfluids see, 
for example, [Q). 

Let us discuss the consequence of the vanishing of the 
uiq moment (^). With assumption ([ij) for the strength 
distribution one immediately finds the result a~^—a^ and 
the splitting between the two frequencies can be directly 
written as 



— 1712 /n^t 



^2 = {[[F^,H], [H,F+]]) = 7V([j_, , (8) Use of dH-O) then yields the relevant result 



where the average ( ) is taken on the state |0) which 
may or may not contain a vortex and we have used the 
property f| = F_ . 

The first commutator (^ vanishes because the opera- 
tors F^ and F_ , depend only on the spatial coordinates. 
The double commutator (^ is the analog of the f-sum 
rule and gets contribution only from the kinetic en- 
ergy term since both the external potential and the two- 
body interaction commute with F±. Finally the current 
operators j± = [p^ /2M, f±\ entering the third sum rule 
are defined by 



2 _ 7LJ^K f Na\ 

M {rl) A2/5 \ ^ a^J 



-2/5 



(13) 



(14) 



for the TO = ±2 modes, where a± = ^Jh/Muj±_ is the os- 
cillator length in the radial direction, while A = Wz/w_l 
characterizes the deformation of the harmonic trap. The 
same calculation can be carried out for the modes excited 
by the m = ±l quadrupole operators (||). In this case the 
result is 



J± 



2Mi 



V/±(r).p + /i. 



(9) 
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(15) 



where p is the usual momentum operator. Evaluation of 
t2 is straightforward in the case 



the sum rules to^ and to 



of the quadrupole operators (||) and 
find the result 
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For TO = ±2 we 

(10) 
(11) 



where Iz is the zth component of the angular momentum 
operator and k is the quantum of circulation of the vortex 
(k = ± 1, ± 2, ...). Of course in the absence of vortices 



In the last equalities of dl^ p^ we have used the Thomas- 
Fermi approximation [|15[ to evaluate the square radii of 
the condensate. Notice that for spherical trapping the 
splitting between the m = ±2 frequencies is twice the 
splitting between the to = ±1 modes. For large N the 
shifts become smaller and smaller showing that that in 
this limit the effects associated with the current of the 
vortex are small corrections to the collective flow of the 
oscillation. Neverthless the splittings can be sizable. For 
example using a spherical configuration with a/a± — 10"'^ 
and = 10^ and taking one quantum of circulation (k = 
1) one finds that the relative shift (w+ — uj-)/uj of the 
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m = ±2 states is about 10%, having used the large N 
result \^L0± for the average frequency oj. This shift is 
much larger than the typical experimental uncertainties 
in the measurements of the collective frequencies 

It is worth pointing out that the frequency shift due 
to the vortex is exhibited by the quadrupole excitations, 
but not by the dipole modes excited by the operators 
f± = X ± iy. In fact in the dipole case the operators 
j± = }i{Pxay)/Mi commute and the sum rule iden- 
tically vanishes. As a consequence one finds a;+= lo^ as 
expected from general arguments. Notice however that 
in addition to the above modes excited by the center of 
mass operators, another dipole mode, localized near the 
core of the vortex, has been predicted to occur with 
frequency different from uj±. It has been recently sug- 
gested [Im that this mode could play an important role 
in driving the instability of the vortex. 

The above results for the shift of the quadrupole fre- 
quencies hold for large N, where the assumption that 
the operators /+ and /_ excite a single mode is justi- 
fied. When the adimensional parameter Na/a± becomes 
small, this assumption is no longer valid. In particu- 
lar, in the limit of a non interacting gas, a vortex with 
quantum circulation k=+1 corresponds to putting all the 
atoms in the 1 p state (Z^ = +1) of the harmonic oscillator 
hamiltonian and the m = — 2 operator /_ can give rise to 
duj — 2LUQ as well as to 610 = 0100 transitions. Viceversa the 
operator /+ gives rise only to dLo = 2LOQ transitions. The 
corresponding strengths are, respectively, = 2a'^N, 



Result 



for 7714 provides a crucial check of the valid- 



•^down = 4ai7V and <j+ 



6a\N. 
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moments m 
the quadrupole operators. For example, for m 
find 



In order to study the transition from the noninteract- 
ing to the large N regime we have to remove the single- 
mode assumption (|l^ and hence we need the knowledge 
of additional moments of the distribution function. The 
and can be easily calculated for 

± 2, we 

(16) 
(17) 
, (18) 
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-£^kinj_ 



M2 

32tfuj 



-N{r' 



l^LO^. 



1 



int 



i6£;ho. 



where i^kin^ and -Eho^ are the radial contributions to the 
kinetic and oscillator energies respectively, and 
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dr 



no 



|V±no| 



no 



no 



is the contribution to arising from two-body interac- 
tions where no is the density of the condensate and k is 
the quantum of circulation of the vortex. Result (|l^) for 
ni^ permits to obtain directly the asymptotic behavior 
of the quadrupole frequency. In fact in the large N limit 
the kinetic energy term in (^6|) is negligible and the ratio 
(m^ /m^y^'^ approaches the hydrodynamic result ^/2io±. 



ity of (|lj). In fact one can see that in the large N limit 
the single-mode approximation (|T^, with the dispersion 
law 



UJ± 



= V2uj± ± Aa 



(19) 



where Alo 
result (|T^) 



= — iO- is given by ( p^ ) , is consistent with 
for the sum rule m^, up to effects linear in 



Alu. Differently from the other sum rules, rrig depends 
explicitly on the two body interaction. This contribution 
is very sensitive to the core region of the vortex and is im- 
portant to describe the crossover from the noninteracting 
to the Thomas-Fermi limit. 

We have calculated numerically the sum rules to^ with 
p= 1, .., 5 using the solution of the Gross-Pitaevskii equa- 
tion for a vortex of quantum circulation k = -1-1 [17j. 
These sum rules are then used to evaluate the quadrupole 
energies and strengths, by assuming that the operator /_ 
excites two states with different energy, in analogy with 
the behavior of the noninteracting model. The results 
for a spherical potential are reported in Figure 1. where, 
for semplicity, we have plotted only the low frequency 
mode excited by /_. The strength relative to the high 
frequency mode excited by /_ in fact vanishes rapidly 
when N increases and hence this mode is not physically 
relevant, except for small values of N. In the Figure we 
also report the average energy {m^ /m^Y^^ calculated in 
the absence of vortices. This energy turns out to be very 
close to the exact numerical solution of the linearized 
Gross-Pitaevskii equation (differences are always smaller 
than 0.5%). Notice that the frequency w^q^^^ corresponds 
to a rigorous upper bound to the lowest frequency of 
modes excited by f±. 

The frequencies lo± start respectively from 2loj_ and 
Olo± and approach the value ^/2io± for large N. In the 
same limit the strengths tend to the asymptotic value 
a+= a- = 4:V2 a\N\^/^(lbNa/ai_f''^ /7. The Figure 
shows that for large Na/aj_ the asymptotic dispersion 
relation (|l^ ) well reproduces the behavior of the two col- 
lective frequencies. Notice that for such values of Na/a± 
the strengths tr^ practically coincide, confirming the va- 
lidity of the single mode assumption (|l2|). The small 
asymmetry of the calculated frequencies with respect to 
the hydrodynamic value is a consequence of the fact that 
for the values of Na/a± reported in the figure the exci- 
tation energy differs from \/2lj±_ even in the absence of 
vortices. 

We have also explored the case of negative scattering 
lengths. In this case we find that the excitation energy 
of the lowest mode becames negative revealing that the 
vortex is unstable against macroscopic fluctuations of the 
density. The dispersion law for small value of the param- 
eter a is found to be 



down 



fNa\ /T 



(20) 



and changes sign with a. 
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It is finally interesting to discuss how the above re- 
sults are modified by changing the geometry of the prob- 
lem. This is important because toroidal configurations 
are expected to suppress the mechanisms of instability 
of the vortex A useful way to pin a vortex might 

be achieved with a thin laser beam stabilizing the core 
of the vortex along the z-axis. The thickness d of the 
beam can be a few microns. This is larger than the size 



of the core, fixed by the coherence length ^ 



but can be significantly smaller than the size of the con- 
densate R. In this case the structure of the core of the 
vortex will be significantly modified by the pinning, but 
the macroscopic behaviour of the collective excitations 
and in particular results ( p^jl^ ) for the splitting, will be 
modified only in a minor way. An estimate of this effect 
can be obtained by calculating the change of (r^) due to 
the presence of the repulsive potential generated by the 
laser beam. We expect small effects ii d R. 

A quite different behaviour is achieved by choosing a 
ring geometry. In the case of an ideal ring of radius R the 
problem is analytically soluble using Bogoliubov theory 
p6| . In this case the natural excitation operators have 
the form f± = exp{±im(f)) where is the azimuthal angle 
and ±m is the angular momentum carried by the exci- 
tation. The kinetic energy operator can be replaced by 
(-?iV2Mi?2)a7a02 ^jjg 1^ g^^^ (|) 

come m'^ = h^m^/MR^ and = 2h^ {I / AP R"^ . 



Using the Bijl-Feynman ansatz (12) one immediately 
finds tj+ — LO- = 2hmK/ M R^ , where k is the quantum of 
circulation of the vortex, in agreement with the results 
recently discussed in . Notice that for large R the low- 
est collective excitations in the ring geometry correspond 
to one-dimensional compression waves with dispersion 



R 



± 



h K \ 



MR^ 



(21) 



where c is the sound velocity. These frequencies, which 
are the analog of (^9|), coincide with the ones calculated 
for a system at rest in a frame rotating with angular 
velocity uj = Kk/MR^. 

After completing this paper we re ceived a preprint b y 
A. A. Svidzinsky and A.L. Fetter ( |cond-mat/980318lD , 
based on a hydrodynamic description of normal modes. 
Similarly to the results of this work predicts a fre- 
quency shift also for the dipole mode. The origin of this 
discrepancy remains to be understood. 

Useful discussions with F. Dalfovo, A.L. Fetter, M. 
Gunn and N. Wilkin-Gunn are acknowledged. This work 
was supported by the BEG advanced research project of 
INFM. 
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FIGURE GAPTION: 

Frequencies (a) and strengths (b) relative to the m — ±2 
quadrupole modes in the presence of a k = 1 vortex, 
as a function of Na/a± for a spherical trap ((j^ho — 
flho = A = 1). The dotted lines correspond to the 
large N behavior (|l^) . The arrow indicates the Thomas- 
Fermi limit u) — \/2uj±. Dashed-dots correspond to the 
ratio (m^ /m^y^'^ without vortex. Strengths are given 
in imits oi a*, N. 
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